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SUMMARY 


The purpose of this study was the examination of an empiri- 
cal method of obtaining aircraft transfer functions directly from 
in-flight records without recourse to lengthy frequency domain 
techniques. 

The proposed method utilizes a precision analog computer 
and plotter to generate response curves matching the in-flight 
records. 

The computer set-up used for simulation is so arranged 
that each numerator coefficient of the transfer function is directly 
related to a potentiometer setting. At the most, two potentiometer 
settings are unknown, since the first and last coefficients can be 
obtained from the initial slope and steady state value by using 
Laplace initial and final value theory. In the case of longitudinal 
motion the period of the phugoid is of such a long duration as to 
make it impractical to obtain final value steady state conditions 
in actual flight records. This is not prohibitive for this problem 
Solution method, however, since one coefficient can be determined 
from the flight record initial slope or suitable accelerometer 
measurement and then the two remaining coefficients are de- 
termined by adjusting potentiometer settings. 


Utilization of this method indicated that errors would be 


ii 





of the order of five per cent in magnitude with the combination 
of computer and recorder used. It is felt that with a high speed 
repetitive computer and a large high-persistence cathode-ray 
oscilloscope, the time of solution and degree of error might be 
reduced if sufficiently accurate flight data were obtained 
initially. 

Further study is indicated for inputs differing from step 
functions such as input wita the Laplace transform: 


P. EM 


5(57+ 1) 


The theory for this case is mentioned briefly in the Proposed 


Technique section. 
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INTRODUCTION 


Considerable theoretical work has been done with frequency 
domain techniques as applied to aircraft dynamic analysis, whereas 
little effort is indicated in the field of transfer function determination 
directly from in-flight transient response data. The purpose of this 
thesis is to investigate a simplified method of determining transfer 
functions by the latter means. 

The empirical transfer function solutions are checked 
against the solutions of the equations of motion utilizing stability 
derivatives and other data of the F-86D aircraft. This aircraft 
was chosen Since it is of relatively recent vintage and has been 
in existence long enough to supply reliable information sufficient 
to verify the techniques developed in this thesis. Transfer 
function poles are readily determined from flight test records; 
however, the zeroes require tedious computation when frequency 
domain methods are used. The values of the transfer function 
poles were taken from analytical computations, however, they 
were checked on the computer for validity and accuracy. The 
previous precaution was taken in order to narrow the field of 
investigation and source of error in determining the zeroes which 
were the crux of the overall problem. It was subsequently de- 
termined that this precaution was not necessary and would not 


be required in an actual practical situation. 





Two basic theorems of the Laplace transformation (Sec 
Ref. 1) constitute an important part of the development: 
(a) Final value theorem - If the function f(t) and its first 


derivative are Laplace transformable and if 


ШЕСІ = F(s) 


then 
Limit. sF(s) = limit f(t) 
5-0 t—» oo 
(b) Initial value theorem - If the function f(t) and its first 
derivative are Laplace transformable and if 
г [50 | = F(s) 
then 
Limit SF(s) limit f(t) 
Ss —>00 t —> 0 


The application of these theorems will be shown in con- 
junction with measurements taken from the assumed in-flight 
transient response curves. 

This study was undertaken by R. W. Ulgen and W. Р. 
Vosscler under the guidance and direction of Professor R. M. 
Howe as a joint thesis project to satisfy, in part, the require- 
ments for a Master of Science Degree in Aeronautical Engineering 
at the University of Michigan, Ann Arbor, Michigan, during the 


period September 1957 - June 1958. 
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SYMBOLS AND NOTATION 


Angle of attack, degrces or radians 

Wing span, fcet or inches 

Chord length, feet or inchcs 

Drag coefficient, D/qS 

Lift coefficient, L/qS 

Momont coefficient about aerodynamic center 
Normal force coefficient, N/qS 


Normal force cocfficient curve slope duc to 
angle of attack dCy/do 


Pitch moment coefficient, m/qSd 

Slope of the lift curve 

Elevator deflection, degrees or radians 
Acceleration of gravity, ft/sec” 
Altitude, feet 

Lift, pounds 

Mass of airplane, slugs 

Mach number 

Dynamic pressure, jo v^ /2 

Wing area, squarc teet 

Airspeed, velocity, knots, mph, ft/sec 
Weight of aircraft, pounds 

Angle of attack in yaw planc, degrecs or radians 
Natural frequency, rad/sec 

Damping coefficient 


Time constant, scc 





AIRCRAFT AXIS SYSTEM 


Axis 

Force 

Moment 

Linear Velocity 
Linear Acceleration 
Angular Displacement 
Angular Velocity 
Angular Acceleration 


Moment of Inertia 


iU*jV* kw 


= iP-jQ-tkR 


Qi QU 
| 


E > 
= +] м 
ZH N 


u V W 
u V W 
ф ө y 
p q r 
р а ; 
2. D I, 


space linear velocity vector 


space angular velocity vector 





ILLUSTRATIONS 
1 Equipment 
2 Longitudinal Computer Diagram 
3 Longitudinal Response Curves 
4 Longitudinal Simulator Diagram 
5 Longitudinal Simulator Response Curves 
6 Lateral Computer Diagram 
7 Lateral Response Curves 
8 Lateral Simulator Diagrams 


3 Lateral Simulator Response Curves 





PROPOSED TECHNIQUE 


THEORY 
Transfer functions can be obtained from the equations of 


motion which are ratios of polynomials of the general form: 


n-1 
хо). bas” + 15 + — bo 96 
f(s) | àmsP + apm-1sSm-l 4 —.......ag 


Multiplying the above by 8171 апа employing the Laplace 


initial value theorem for a step input: 


m-ny b 
ету е А пс 
S — 00 t —»0 dt am 


The slope of the (m-n-1) derivative response curve is then 
equal to b, б/а: 
The final value theorem for a step input gives: 


limit sX = limit X(t) = 2o 
S —»0 [> од О 


The first and last coefficients of the numerator can now 
be easily evaluated since the denominator coefficients are found 
by the simple means of measuring damping and frequency. For 
n 2 2 there are (n-1) unknown coefficients remaining. 

If the above transfer is mechanized (n-1) potentiometers 
representing the unknown coefficients can be varied so as to 
obtain results matching the in-flight records. This is not overly 
difficult to accomplish since there are a maximum of two 


potentiometer settings which are unknown and these are known 





at least to an order of magnitude. This would apply to rigid- 
body aircraft. 

A suggested method of procedure for a step function in- 
put is to match phase fairly closely and then match initial slopes 
for the same Steady state value. Iteration brings relatively 
rapid convergence. 

For the case where the input transform d (s) is 4 

s(s 1) 
the final value theorem gives results identical to the previous case. 


The initial value theorem, however, indicates that 


limi ехо 209 


t 0 dem-n-1 TY am 
where 7 is the time constant of the input. 


Only examples having step function inputs will be con- 


sidered in this thesis. 





EQUIPMENT 

The computer used for these studies was a 20 amplificr 
electronic differential analyzer designed and built by the De- 
partment of Aeronautical Engineering, University of Michigan. 
A photograph of thc computer, along with the Variplotter, is 
shown in Fig. 1. This computer had resistor accuracy of . 02%, 
which was more than required for the present study. The 
computing amplificrs had a gain of approximately 10° and were 
manually balanced. Coefficicnt potcntiometers were sct to 
an accuracy of 0.1%. 

The Pace Associates Model 1100D Variplotter used to 
plot the analogue computer outputs is a tablc-top size, self- 
balancing, potentiometer type of recorder. It records two 
variables simultaneously in thc form of a rectangular coordinate 
graph. The static accuracy of both the pen and arm is . 075 per 
cent. The dynamic accuracy for both is . 2% for plotting speeds 


up to 15 inches per second. 
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EXARMELCES 

Examples of the determination of zeroes for both longi- 
tudinal and lateral functions will be given. 
I. Longitudinal Transfer Function Zeroes 

Since no suitable flight records were available, flight 
data was assumed from solving the aircraft equations on an 
electronic differential analyzer. 

Under conditions of perturbations from steady level 
flight at Mach .8 and 30, 000 feet the following equations 


were utilized based on the airplane data in Appendix A. 


ù = -.00664u + .0497w - 31.9, - 32.20 

W z -.0407u - 1,018“ + 790.34 - 1.29 0 - 1.55 B. 
q = .000475u - .0119w - 137 q - .413 б. 
ө = а 


Scaling the set of equations for use on the analog computer 


gives 

ü = -.00664u + .00497(10w) - .319(100q) - .322(100 O) 

W = -.0407u - .1018(10w)+7.903(100q) - .0129(1008) -. 155(1003 
1004 = .0476u - .119(10w) -.737(100q) -. 0413(1092) 

100 - 1004 


These equations were then mechanized as shown in Fig. 2 
(see Appendix B for all remaining figures). The following scaling 
was used: 
l fps = .1 volts 


l rad = .1 volts 





1! 


1 deg EP volts 
lrad/sec = .1 volts 
The exact transfer function of Yes was found from 
e 


this set of equations to be: 


2 
В + .983 5 - 00855 


+ 1.76253 + 10.165“ + .01958 +. 0311) 


Since there is very little coupling between the short period and 
the phugoid motion it was considered possible that the poles and 
zeroes of the above transfer would be the same as the poles and 
zeroes of the uncoupled transfer functions. 

The approximate formula for the short period transfer 


function as based on Ref. 2 is 


Kg fe (1.032 s + 1) 
5 (. 0978 52 + .1885 5 +1) 


44 


е 
The approximate formula for the phugoid transfer function 

as obtained from the same source is 

KƏ Se (а А 1) 


9/4 м 
© (305.5 52 + 2.4058 + 1) 


The product of the above two expressions does indeed give 
a very close approximation to Y e Ф, which makes a lengthy 
Newton Raphson approach (see Ref. 3) or other factorization 
technique unnecessary. 

The parameters of the approximate transfer function 


are then as follows: 
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,0572 radians 


= 3.22 radians 
а) ns Sup 


s = 1.862 seconds Тр 


ie = .304 УЬ 


These values were used to assure that thc rcsponse 


T 


109.5 scconds 


. 0685 


curves obtaincd from the mechanization of these curves were 
indecd correct. The computer rcsponse curves of Fig. 3 had 
the following parameters 


(9 3.22 radians (А) = ,0572 radians 


ns 


i 
5 
O 


T 2 seconds T 112 seconds 


3 7 p 


£ 


The above paramctcrs give the poles immcdiately in the 


20 


| 
co 
TUS 
= 
ї 


Canonical form. 
The determination of zcroes for the transfer function 
Yo was considered next. 
бе 
In the interest of a simpler circuit, the constant term 
(.0311) in thc denominator was assumed to be negligible. An 


's' was then factored out of the remaining terms giving the transfer 


function of Y 
ade 


2 
-23.6(s? + .983s + .00855 
МИА = A чо)... 


е (5% + 1.76252 + 10.168 + .0195) 
which has the dimensions of degree per degrce or radian per 
radian. 
The denominator is now a cubic with coefficients which 
can be equated to potentiomcter scttings and amplificr gains as 


follows for the longitudinal simulator computer circuit of Fig. 4: 
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“кук, = 1.762, Ki EUIS Ко - "EU; 1 = · 1762 
Ку Кока = 10. 16, Ko = 2) p = „008 
ОК Кокакг = e 0915, (кк). 22 . I; Q3 == „ОГО 


Inasmuch as the steady state of the phugoid is difficult 
to obtain in flight tests, it was decided that a more practical 
method of procedure would be to determine the potentiometer 
setting affecting the initial slope by utilizing flight records and 
to vary the remaining two potentiometers. 

The initial slope potentiometer setting corresponding to 
the coefficient of the s^ term in the numerator was found in the 
following manner. 

The initial slope is found from Fig. 3a or 3b. Recalling 
that 

A Iure bo/ an 
we have from Fig. 3 
Бо = .4 rad/sec“/degree. 

This was used to obtain the potentiometer setting as 
shown in Appendix C. 

As has been indicated before, steady state for the phugoid 
takes too long to obtain and a cut and try approach is more 
expedient. This leaves potentiometers five and six which must 
be adjusted in the computer setup of Fig. 4. Various combinations 
over the full range of each potentiometer were tried. At each 


setting an output trace of q for a 1° elevator step was compared 
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with the assumed flight record. Noting the effects and trends of 
these combined settings it was possible to come within 5% of 
what the actual required settings should be after only about 18 
trials. Increased magnification using specialized high-persistence 
cathode ray oscilloscopes should increase aecuracy of sueh 
measurements. 

Fig. 5 illustrates the technique in a general way uti- 


lizing the following sequenee of settings: 


Trial Potentiometer 
9 6 
А . 5 5 
В 219 9 
С 1::0 „9 
р 1.0 ‚0 


It is to be noted that the initial slope remains unaffected 
as might be expected. The latter condition suggests that the 
potentiometer 4 setting could also have been obtained by trial 
and error matching instead of from initial slope measurement. 
Variation of potentiometer 5 affects the time and magnitude of 
the first peak as well as setting the continued level of the entire 
eurve. Variation of potentiometer 6 has a negligible effect on 
the initial peak, however, it raises or lowers the tail of the 


curve thus influencing the final or steady state value. 





II. Lateral Transfer Function Zeroes 
Again, since no suitable flight records were available, 
the flight response curves were obtained by solving the air- 
craft equations on an electronic differential analyzer. 
In a manner similar to the longitudinal case, the air- 
plane data in Appendix A was used to obtain 
У = -. 179v + 31.8p - 796г + 32.2 Ф 
р = -.435v - 3.6p+ 1.295r + .313r - 1.015 Eh 1644. 
t= .0121v - . 0389p - . 388r + .0765p - . 0146 б - „0838 > 
2 -p+.04r 
Scaling the set of equations for use on the analog 
computer gives: 
v--.179v +.318(100p) - 7.96(100r) + . 322(100) 
100р --.0435у- 3. 6(100р) + 1.295(100г) + .313(100г) - 101.54. 
+ 16. 4d, 
100r = 1. 21v - .0389(100p) - .388(100r) + .0765(100p) - 1. 46 d, 
-8.38 $, 


1002 - 100p + .04(100r) 


These equations were then mechanized as shown in Fig. 6. 
The scaling was as follows: 

l fps = . 1 volt 

l rad = . 1 volt 


1 deg = .1 volt 
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A static check was performed by observing the voltage 


on inverters 1, 6, and 11 for various initial condition voltages. 


The response curves are shown in Fig. 7. 


The important transfer functions for this set of equations 


were found to be as follows: 


ty >. 2 - 
Уф Or = .138(s^ - 2.455 - 9.86) 


D(s) 


2 
YD Sa - 1.05 0825 + 10.1) 


D(s) 


M - -62.7(82 + 4. 068 – .0296) 
vor 


D(s) 


-42(s^ - .872s -.325) 


Y 
vba 5.2 


Y = .073(s?-- 4.4982 4- 1.4568 + .752) 
rdr 


D(s) 


E 3 2 
Y da - .0945(s" t .321s" + 4.48s + 4.51) 


D(s) 


п 


where 


D(s) = s* + 4.1655 + 10.652 + 38.58 - . 01975. 


Since the transfer functions for uncoupled rolling motion, 
Spiral motion, and the dutch roll motion promised to be useful in 
finding the poles and zeroes of the above transfer functions these 


modes of motion were investigated. For uncoupled rolling motion 


we assumed that r = (5 = oe = 0. This gives 


р = -3.6р - 1.015 df, 
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or 

ф o . 382 | 

Sa 8(. 2878 + 1) 

For the dutch roll motion we assumed that p - ф = 9, - 0. 
This gives 

v = -.179v - 796r 

г = .0121v - .388r - .0838 d, 
Or 


MS, = 24 28075 89:72) 

For spiral motion we assumed thatv = p = r = 3 - d, - 
d. - 0. This gives 

0- -.179v + 31.8p- 796г + 32. 2 ф 

0 = -4.35у - 360p T 129.5r 


0 = 1.21v - 3.89p - 38.8r 


О ф -p- .04r 
If we assume in addition thatv= p = 0, we get a spiral mode 
time constant of 1875 seconds. 

With these as approximations, synthetic division coupled 
with the Newton-Raphson iteration technique for determining roots 


was used to give 


Үф4. 


-.138(s - 3.267) (s + 3.022) 
D(s) 


2 
Үф «Фа = 10518 + .5828 + 10.1) 
D(s) 
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“ба. T(8"-- 4. 06) s 
D(s) 


ње 

“ 

zn 
H 


x E = cdd SEI пне) 
уба 


D(s) 
ү р = .073(5 + 4.1) (52 + .32s + .2) 
ГГ 
D(S) 
Y d. - .0949(s + .91) (s - .579s + 5) 
га 
D(s) 


where 


Г(в) - (5 - .000513) (в + 3.95) (s + . 22555 + 9.742) 


The parameters associated with the denominator of the 


various uncoupled transfer functions are then as follows: 


Тр = :287 seconds - the time constant of the uncoupled rolling 
motion. 
T = 1875 seconds - the time constant of the spiral mode 


w оар = .324 radians per second - dutch roll 
Mn - .875 

The parameters associated with the denominator of the 
transfer function are found to be: 

цэ = „2598 seconds 


T. = 1950 seconds 


S 
IDE 
Par = ‚875 


These parameters give the poles immediately in the 
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canonical form. 
The parameters of the simulator circuit shown in Fig. 8 
are determined by setting the coefficients equal to the potentiometer 


settings and amplifier gains as follows: 


x, KK, - 4,16 Ky - 1, Ко - 10 21 = .416 
2 = ex - 

К.К, = 1006 Ко = 2, Ээ. 
Cx a = = 

ак Кокак, = 98.4 Ko = 2, cx цг 96 


since 


D(s) = s(8° + 4.1682 


+ 10.68 + 38.4) 

Two transfer functions, Yos. and Уг fr were investigated, 
as being representative. The value of the potentiometer setting 
corresponding to the highest degree 's' term in the numerator was 
found as in the longitudinal case. See Appendix C for details. The 


remaining potentiometers were adjusted giving the plots and variations 


Shown in Fig. 9. 





20 


RESULTS AND CONCLUSIONS 

By determining potentiometer settings corrcsponding to 
initial slope and also the steady state valuc, where available, 
it was possible to obtain accuracy in matching the curves such 
that the zeroes were determined with errors of lcss than 5%. 

In the case of q/f, , where two potentiometer settings 
were unknown, only 18 trials were necessary to obtain this 
accuracy. 

See Fig. 5 and Fig. 9 to observe the variation in re- 
sponse curves for variation in numerator potentiometer settings. 
It is concluded that the best general procedure is as follows: 

1) Obtain estimates of poles and zeroes from pre-flight 

data. 

2) Obtain poles from flight records involving step inputs 
and coefficients of first and last terms in numerator 
corresponding to initial slope and Steady State, if 
possible. Then compare with previous estimates. 


Limit amarga P = b Mo Иа 
t—» 0 


Limit X(t)- bó /a, 
t —= 00 


3) Adjust potentiometers to correspond to the known 
parameters. 
4) Set estimates in those potentiometers which must 


be adjusted. 





9) 


6) 


8) 


9) 


Make preliminary adjustments to observe which 
potentiometer has more effect in the case where 
two potentiometer Settings are unknown, 

Adjust the more critical to give loose matching, 
observing trend of frequency, phase shift, peak, 
and steady state values. 

Adjust the remaining potentiometer for a closer 
match while observing trends. 

Return to the first potentiometer and repeat ad- 
justment. 


Continue iterating for as close a match as desired. 
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APPENDIX A 


AIRCRAFT DATA AND STABILITY DERIVATIVES 


The following data is given in references 4 and 5 for 


the F-86D aircraft at 30,000 feet, Mach 0.8, weight 15, 000 lbs. 


S = 288 ft? 
I = 8625 slug ft“ 
XX 
2 
12-15 300 slug ft 
V = 796 ft/sec 
p 
& = 2.32 20.04 rad 
3 
/2 = 00089 slug/ft 
Сү, - 0.185 
2156.5 
m 6 
Cra - 4.62/radian 
C 20! 
Dc, 


E = -0. 424 
Мы 


Ом, = -1.5 
Сре = „ 143 
Ср = . 018 


MD - -.0185/degree 


С = 0. 0089/degree 
aie = 


с = 8.09 ft b = 37.12 ft 
lyy = 29400 slug ft 
I, = 2, 700 slug ft^ 
Ee = -,82/radian 
ыг = 
E = 
шү men 
Cy, za 
8 
Сү, --9.10/гайап 
“Ng = 117 
Cp, = =. 433 
S =. №9 
CNp s cU 
Exa E 
E - -.0029/degree 


C = .00043/degree 
ID / deg 


Cy 5 = - 000055/degree 


CN 5, = = 001 


“Ng sally 


C = -, 10 
L 
/8 











| 

| o 

| = = 
| 





APPENDIX B 








-lOw 


(ООа 


1008 


ІО м 
-1004 
1008 


| Обе 


lOw 


100 q 


О Фе 


(v) (=) 8 
О о 
О бу 
o > 
7 8 

СЛ 


2) ; 
ni 
ce 
| 
= 


Со 


Ж 
> 


O 
N 


JOI8 I 


19 1 


23 


(5) zw -lOw 
Qu | 
| | | 
0476 | 
190 | 
| | 
737 | й 
È 100 q -1004 ооа 
413 | 
© | 
: | 
-100q=-1008 ioo 8 
Fig. 2 


Cumpufer Circuit for Solving the Longitudinal Flight Equations 
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APPENDIX. 








APPENDIX C 


Determination of Potentiometer Settings 
For the longitudinal case we found that 
bo = .4 rad/sec^/degree. 
from Fig. 3 a or 3b. 
The scale that we wish the plotter to have is 
1 inch = .1 rad/sec“ 
with the pen scale factor set to 10 volts per inch. 
For a one degree step input the initial q should be equal 
to 40 volts, since 
40 volts = .4 rad/sec”. 
The equation relating the 40 volts to the parameters of 


the circuit can be determined by reference to the figure below. 
2. 
Ves 


95 озбауы) 





«——— [2/»ts)] (4K, 5e(s)/kgk7) о = ХО 
Limit sX(s) = 4, = (2724) («КБ d. / KK) - 40 volts 
8--» ОО 


For d Q7 2 volts 
«К - 2 
In brief then for the initial slope 
(К/а,) (ё6,К, 9 ¿/KgK7) = (5, ¿/K i) (09/89) 


and for steady state 


(K/a,) (Ki Dg) = (Kyi/Ky;) (6/20) 





K = denominator multiplying factor 


K = the pen scale factor on the plotter (in the above 


example it was ten volts per inch. ) 
ESA - the number of radians /second or the number of 
radians/second /second (in the above example 


it was . 1) 
































